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Abstract 

This article is devoted to the regular fractional Sturm-Liouville eigenvalue problem. Ap- 
plying methods of fractional variational analysis we prove existence of countable set of orthog- 
^ ■ onal solutions and corresponding eigenvalues. Moreover, we formulate two results showing 

OO ' that the lowest eigenvalue is the minimum value for a certain variational functional. 
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Introduction 



In 1836-1837 French mathematicians Sturm (1803-1853) and Liouville (1809-1855) published series 

C^ , of articles initiating the new subtopic of mathematical analysis - the Sturm-Liouville theory. It 

deals with the general linear, homogeneous second-order ordinary differential equation of the form 

^ (p(a;)^j +9(a;)2/ = Aw(a;)y, (1) 

where x € [a, 6] and in any particular problem functions p{x), q{x) and w{x) are known. In 
addition, certain boundary conditions are attached to equation ([T]). For specific choices of the 
boundary conditions, nontrivial solutions of ([T]) exist only for particular values of the parameter 
A = A'"'\ m = 1,2,.... Constants A*^™^ are called eigenvalues and corresponding solutions y(™)(a;) 
are called eigenfunctions. For a deeper discussion of the classical Sturm-Liouville theory we refer 
the reader to [51[ini[12- 

Recently many researchers have focused their attention on certain generalization of Sturm-Liouville 
problem. Namely, they are interested in equations of the type ([1]), however with fractional differ- 
ential operators (see e.g., pQ [H [HI UHl Ell [HI [131 [13 [ISl [H]) which are integrals and derivatives 
of arbitrary real or complex order. As mentioned in [3], Science Watch of Thomson Reuters 
identified the fractional calculus as an Emerging Research Front area. This is due to its many 
applications in science and engineering. A comprehensive study of the fractional calculus and 
its applications can be found in several recent books [H [ll [3 [l7] . In this note we develop the 



fractional Sturm-Liouvillc theory by studying the Sturm-Liouvillc eigenvalue problem with Ca- 
puto fractional derivatives. We shall show that fractional variational principles are useful for the 
approximation of eigenvalues and eigenfunctions. Traditional Sturm-Liouvillc theory does not 
depend upon the calculus of variations, but stems from the theory of ordinary linear differential 
equations. However the Sturm-Liouvillc eigenvalue problem is readily formulated as a constrained 
variational principle, and this formulation can be used to approximate the solutions. We empha- 
size that it has a special importance for the fractional Sturm-Liouville equation since fractional 
operators are nonlocal and it can be extremely challenging to find analytical solutions to fractional 
differential equations. Besides allowing convenient approximations many general properties of the 
eigenvalues can be derived using the variational principle. 

The paper is organized as follows. In Section [T] a brief review of the fractional calculus is given and 
three technical lemmas are shown. Our main results are then formulated and proved in Section [2 
we show the existence of orthogonal solutions to the fractional Sturm-Liouville eigenvalue prob- 
lem (Theorem[5]), prove that lowest eigenvalues can be characterized as minimum values of certain 
functionals (Theorem [3] and Theorem|4]). In Subsection l2.3l we illustrate our results through an ex- 
ample. Finally in Appendix, reader can find two lemmas concerning certain convergence properties 
of fractional and classical derivatives, that play an important role in the proof of Theorem [21 

1 Preliminaries 

The reference books for the fractional calculus are [6l [71 [H [17] . Here we only recall necessary 
definitions and properties of fractional operators. Moreover, we prove lemmas that will be used 
in the proof of the main result. Theorem [2 

Definition 1 (Left and right Riemann-Liouville fractional integrals). We define the left and the 
right Riemann-Liouville fractional integrals i^^ and J^_ of order a CzR (a > 0) by 

and 



r(a) Jx (t - x) 

respectively. Here T{a) denotes the Euler's gamma function. 



^-/(-) ■■= Wit / u'!Z. ^ ^ e [a, 6), (3) 



The following assertion shows that Riemann-Liouville fractional integrals satisfy semigroup prop- 
erty. 

Property 1 (cf. Lemma 2.3 0). Let q;,/3 > and f G LP{a,b), (1 < p < ooj. Then, equations 



i?-/L/(^) = C/w 



are satisfied. 



Definition 2 (Left and right Riemann-Liouville fractional derivatives). The left Riemann-Liouville 
fractional derivative of order a G M (0 < a < I) of a function f , denoted by D"_^_f, is defined by 

yxe{a,b], D:+fix):^Dll:i^fix). 

Similarly, the right Riemann-Liouville fractional derivative of order a of a function f , denoted by 
D^_f, is defined by 

Vxe[a,6), D^_f{x):=-Dllz''fix), 

where D — -r-. 

ax 



As one can see below Riemann-Liouvillc fractional integral and differential operators of power 
functions return power functions. 

Property 2 (cf. Property 2.1 [6]). Now, let a,/3>0, then the following identities hold: 

P^^it-af-\x)^-^^{x~a)^+--\ 
T(p + a) 

D:4t~af-\x)^^^^^^{x^a)^—\ 

i^_ib-tr~\.)^-^i^{b~..r+--\ 

i[p + a) 

and 

D^(6-t)^-i(.x) = -^(6-.r"-i. 
r(/3 - a) 

Definition 3 (Left and right Caputo fractional derivatives). The left and the right Caputo frac- 
tional derivatives of order a G M (0 < a < 1) are given by 

yxe{a,b], ^D:+f{x):^D':+[f{x)~f{a)], (4) 

V.TeK6), ^D^_f{x):=D^_[f{x)-f{b)]. (5) 

Let < a < 1 and / G AC[a, b], then Caputo fractional derivatives satisfy the following relations: 

''D^^fix) = ll+'^Dfix) 

and 

'D^_f{x) = -llz"Df{x), 

respectively. 

Property 3 (cf. Lemma 2.4 0). If a > and f E LP{a,b), (I < p < ooj, then the following is 
true: 

D^+Pa+fix) = fix), 

D^_]^_f{x) = fix), 
for almost all x G [a, b]. If function f is continuous, then composition rules hold for all x G [a, b]. 

The above property shows that the Riemann-Liouville derivative is the left inverse of the Riemann- 
Liouville integral, but we cannot claim that it is the right inverse. More precisely, for 1 > a > 
we have the following situation. 

Property 4 (cf. Lemma 2.5 and Lemma 2.6, 0). // / £ L^ia,b) and ll+"f,llz"f e AC[a,b], 
then the following is true: 

I^+D^^fix) = fix) - t_^il-^(a), (6) 

J^_D^/(x) = fix) I^^^jiz"(6). (7) 

Next results show that for certain classes of functions Caputo fractional derivatives are inverse 
operators of Riemann-Liouville fractional integrals. 

Property 5 (cf. Lemma 2.21 [B]). Let a > and a £ N or a ^ N. If f is continuous on the 
interval [a, b], then 

^D^^r^^fix) = fix), 

'D^^I^_fix) - fix). 



Property 6 (cf. Lemma 2.22 [5]). Letl>a>0. If f e AC[a,b], then 

Pa+'D^+fix) = fix) - f{a), (8) 

i?_^Z?^/(.T)=/(.T)-/(6). (9) 

Note that, if /(a) = 0, then we have 

r,+'D:+f{x) = Il+Dfix) = fix) - f{a) = fix), (10) 

and similarly for /(6) — we obtain 

r.^'Dt^fix) = -ll^Dfix) = fix) - fib) = fix). (11) 

For p-Lebesgue integrable functions Riemann-Liouville fractional integrals and derivatives satisfy 
the following composition properties. 

Property 7 (cf. Property 2.2 [3]). Let a > ;3 > and f G L^ia, b), (1 < p < ooj. Then, relations 

D^,^p,^fix) = r:^^fix), 

and 

DlT^_fix)^l^:^fix) 

are satisfied for almost all x G [a, b] (in the case of continuous f they are satisfied for all x G [a, b]). 
In particular, when /? = A: G N and a > k, then 

D'^r:^fix) = r:-'fix), 

and 

D^r,_fix) = (-1)'= r.z'^fix). 

In classical calculus, integration by parts formula relates the integral of a product of functions to 
the integral of their derivative and antiderivative. As we can see below, this formula works also 
for fractional derivatives, however it changes the type of differentiation: left Riemann-Liouville 
fractional derivatives are transformed to right Caputo fractional derivatives. 

Property 8 (cf. Lemma 2.19 [8]). Assume that < a < 1, / G AC[a,b] and g G X''(a, 6), 
A < P < ooj- Then, the following integration by parts formula holds 



fix)DZ+gix) dx = / gix)^D^_fix) dx + fix)l'-^'^gix)l^^ . (12) 

J a 

Finally, let us recall the following property yielding boundedness of the Riemann-Liouville frac- 
tional integral in the ^^(a, 6) space (cf. Lemma 2.1 0). 

Property 9. Let /3 G M+ and p>l. The fractional integral operator /f ^ is bounded in L^ia, b): 



|/f+/IU.<A>||/|Up K,^^-—^. (13) 



Now, we are in position to prove lemmas that will be used in the proof of Theorem [21 
Lemma 1. Let a G (0, 1) and function 7 G C[a, b]. If 

l-b 

-fix)D''D2+hix)dx = 

for each h G C^[a, b] such that D'^D'^^h G C[a, b] and fulfilling boundary conditions 

hia) = Il+^hib) = 0, (14) 

W^^hix)U^a = 'D^^hix)U^b = 0, (15) 

then jix) = Co + cix, where co,ci are some real constants. 



Proof. Let us define function h as follows 

h{x) := /!+" (7(x) - CO - cix) (16) 

with constants fixed by the conditions 

Il+ hi^) - Co - cix) U^fc = (17) 

4'+ (7(2;) - Co - cix) U^fc = 0. (18) 

Observe that function h is continuous and fulfills the boundary conditions 

h{a) = lll°'h{x)\^^h = /^^ (7(2;) - Co - Cix) \x=b = 



and 



In addition 



We also have 



"^D^^h{x)U=a = D:^h{x)U^a = i?/f+ ilix) - CO - Cix) U=. 
= II+ (7(2;) - CO - Cix) \x=a = 

=D^+/i(x)U=fc = D2+h{x)U=i, = Dll^ (7(2;) - Co - cix) U=b : 
= /q+ (7(2;) - Co - cix) \x=b = 0. 

h'{x) = Ia+{l{x) - Co - cix) G C[a, 6] 
D'^D'^^h = 7(x) - Co - cix e C[a,b]. 



b 

(7(0;) - Co - cix) D'D'^^h{x)dx 



I [-co-cix)D''D2+h{x)di 

J a 



-co • ^DZ+h{x)\%zl. - c,x ■ ^DZ+h{x)\%z''a + ci • 4^T"/i(2;)|^=^ = 0. 
On the other hand 

D^DZ+h{x) = D^DZ+llX^ (7(.2;) - CO - cix) = 7(2:) - co - c^x 

and 



pO no 

0= {■j{x) -co-cix)D''D2+h{x)dx^ (jix) - cq - cix) 

J a J a 

Thus function 7 is 

7(x) = Co + cix. 

Lemma 2. Let a e {^,l) , j £ C[a,b] and D^^"7 e L^[a, 6]. // 

-i{x)D''D2+h{x)dx = 



da;. 



D 



for each h £ C^[a,h\ such that h" G i^[a, 6] and D'^D^_^_h G C[a, 6] fulfilling boundary conditions 
(|14p . ()15p . i/ien 7(.t) = co + Cix, where cg, Ci are some rea/ constants. 



Proof. Wc define funetion h as in tlie proof of Lemma [T] 

hix) := /!+" (7(2;) -co- c,x) (19) 

witli constants fixed by the conditions ([T7]) and (|18p The proof of the lemma is analogous to that 
of Lemma [TJ In addition for the second order derivative we have 

h"{x) = DI"^{-f{x) -co-cix) ^ 

= Dl^" {^{x) -co- cix) = 

Dtr^i^) - (CO + ci«) ^'';°r' - ci 

Let us observe that for a > 1/2 



+ 7 2^ - CO + c.a)'-—-^ c / ; 

1 (a) 1 (a + 1) 



(x — a) 



a-l 



r(«) 

(s - a)" 



eL^[a,6] 



r(a + l) €C[a,fc]cLla,6]. 

Thus, we conclude that h" G L'^[a,b] and function /i constructed in this proof fulfills all the 
assumptions of Lemma [21 The remaining part of proof is analogous to that for Lemma [T] D 

Lemma 3. (1) Let a <E (^, l), functions 7^ e C[a, b], j — 1, 2, 3 and D]^^i, G i^[a, 6]. // 

/ (7i(x)Ma:)+72(xri?:+Ma;)+73(x)I?=i?^+/i(x))dx = (20) 

•/a 

/or each h e C^[a, 6], such that h" G i^[a, b] and D'^D" _^_h G C[a, b], fulfilling boundary conditions 

dUl), HID, thenj3eC^[a,b]. 

(2) Let a G (|, l) and functions 71,2 G (^[a, 5]. // 

{li{x)h{x) + 72(x)^I?^+/i(.t)) dx = (21) 

for each h G C^[a, b], such that h" G L^[a, 6] and D'^D'^ ^h G C[a, fe], fulfilling boundary conditions 
dHI, (USD f/ien 

7i(a:)+'=i^,"_72(x)=0. 

Proof. Observe that integral (|20p can be rewritten as follows 

(7i(x)/i(x) + 72(x-)=7^^+/i(x) + ^3ix)D^D^^hix)) dx = 

(-4V4"-7i(a;) - 4V72(^) + 73(x)) D'D^^h{x)dx ^ 0. 
due to the fact that relations 



i:+ilD^D:^hix) = -h{x) 



and 



li_D^DZ+h{x) = -^DZ+Hx) 
are valid because function h fulfills boundary conditions p4p . p^ . Denote 

7(x) := -Il+I^_^,{x) - Il+i2{x) + l^{x). 



It is clear that 7 e C[a,b] and -Dq+^7 G L'^[a,b]. Thus, accordmg to Lemma [21 there exist 
constants cg and ci such that 

-ll+Ib-ll{x) - ll+l2{x) + -fsix) = Co + CiX. 

Let us note that function 73 is 

73(2;) = ll+I"_-/i{x) + 4^72(2;) + CO + cix. 

Hence its first order derivative is continuous in [a, b] and 73 G C^[a, b]. 
The proof of part (2) is similar. We write integral ([2T|) as follows: 

{ji{x)h{x) + ^2{xrDZ+h{x)) dx = / (-/i+4"_7i(2:) - Il+l2{x)) D^DZ+h{x)dx = 0. 

J a 

The function in brackets is continuous in [a, b] and 

^a;" (-4V4"-7i - 4V72) = -4+4"-7i - 4\72 e C[a, 6] C L2[a, b] 
so we again can apply Lemma [5] and obtain that there exist constants cq and ci such that 

Il+lb-ll{^) + Il+l2{x) = Co + CiX. 

Thus functions 71^2 fulfill equation: 

^i?^ 72 (a;) +71(2^) = 0. 

D 

The crucial idea in the proof of our main result is to apply direct variational methods to the 
fractional Sturm-Liouvillc equation. Starting from the fractional Sturm-Liouvillc equation the 
approach is to find an associated functional and to use this to find approximations to the station- 
ary functions, which are necessarily solutions to the original equation. In the case of the fractional 
Sturm-Liouville equation an associated variational problem is the fractional isoperimetric problem 
which is defined in the following way: 



minJ(2/)= / F{x,y{x),'DZ+y{x))dx, 

J a 

subject to the boundary conditions 

y{a) = 0, y{b) = 
and the isoperimetric constraint 



I{y)= / G{x,y{xVDZ+y(x))dx = ^, 



where ^ G M is given, and 



F : 
G 



[a, b] X K2 
{x,y,u) 

[a, b] X M? 
{x,y,u) 



F{x,y,u), 



G{x,y,u) 



(22) 
(23) 

(24) 



are functions of C^ class, such that S^, ^ have continuous D? derivatives. 

Theorem 1 (cf. Theorem 3.3 [3]). Ify £ C[a,b] with '^D^^y G C[a,b] is a stationary function for 
problem (|22p - (|24p . then there exists a real constant X such that, for H = F + AG, the following 
equation: 



dH 

dy 



ix,yix),^D2+y{x)) + D^^ 



dH 

du 



{x,y{x),^DZ+y{x)) 



holds, provided that 



— ix,y{x),^D2_,y{x))+D^_ 



8G 

— {a:,y{x),^D^^y{x)) 



0, 



^0. 



(25) 



2 Main results 

For the first works on the fractional calculus of variations, we should look back to 1996, when 
Riewe used non-integer order derivatives in order to better describe nonconservative systems in 
mechanics jl9) . Since then, many papers have been written on the subject, for a state of the 
art on the fractional calculus of variations we refer the reader to the recent book [14] . In this 
section, we will present an interesting application of the fractional variational calculus. Namely, 
using the fact that the fractional Sturm-Liouville eigenvalue problem can be remodeled as the 
fractional isoperimetric variational problem, we will show that there exist an increasing sequence 
of eigenvalues and a corresponding sequence of eigcnfunctions for which the fractional Sturm- 
Liouville equation is satisfied. Moreover, we will prove that the lowest eigenvalue is the minimum 
value for a certain variational functional. 

2.1 Existence of Discrete Spectrum for Fractional Sturm-Liouville Prob- 
lem 

We will show that similar to the classical case for the fractional Sturm-Liouville problem there exist 

an infinite monotonic increasing sequence of eigenvalues. Moreover, that apart from multiplicative 

factors to each eigenvalue there corresponds precisely one eigenfunction. 

We shall use the following assumptions. 

(HI) Let ^ < a < 1 and p, q, Wa be given functions such that: p is of C^ class and p{x) > 0; q, Wa 

are continuous, Wa{x) > and (y^Wa)' is Holderian of order /3 < ct — 2- 

Consider the fractional differential equation 

[^D^_p{xrD2^ + q{x)] y{x) = \wo.[x)y{x), (26) 

that will be called the fractional Sturm-Liouville equation, subject to the boundary conditions 

y{a) = y{b) = 0. (27) 



Theorem 2. Under assumptions (HI), the fractional Sturm-Liouville Problem (FSLP) (j26[) - 
(|27p has an infinite increasing sequence of eigenvalues A^^A"',..., and to each eigenvalue A*-"' 
there corresponds an eigenfunction y^'"'' which is unique up to a constant factor. Furthermore, 
eigcnfunctions y'"' form an orthogonal set of solutions. 

Proof. The proof is similar in spirit to [5] and will be divided into 6 steps. As in [5] at the same 
time we shall derive a method for approximating the eigenvalues and eigcnfunctions. 
Step 1. We shall consider problem of minimizing the functional 

J{y) = f [pix){^D:_,yy + qix)y'] dx (28) 

J a 

subject to an isoperimetric constraint 

Hy) = / w^ix)y'dx = 1 (29) 



and boundary conditions ([771) . Fii'st, let us point out that functional (pS)) is bounded from below. 
Indeed, as p{x) > we have 



Jiy)^ I [p{x)CD:^yf + q{x)y'] dx 



q{ Xi f ' q{ x] 

> min -— • / Wa{x)y'^ dx = min — "—— =: A/q > — oo. 

x<^[a,b] Wa{x) J^ x£la,b] Wa{x) 



From now on, for simplicity, we assume that a ~ and b = ir. According to the Ritz method, 
we approximate solution of (j27P " (|29p using the following trigonometric function with coefficient 
depending on Wa- 



(x) = ——= V" Pk sm{kx). 



(30) 



fc=i 



Observe that ym{0) = ym{T^) = 0. Substituting (|5D)l into (pS)) and (P5)l we obtain the problem of 
minimizing the function 

J(/3i,...,/3™) = J([/3])= (31) 



subject to the condition 



, s.cnQ sin(fca;) ^ sin(jx)\ q(x) \ • /• \ 

p{x) "D^^ ^ ' ■ "£'g+ )_' ) + -^ sm(fca;) sm(jx) 



/■iWo 



WQ(a;) 



dx 



/(/3i,...,/3„) = /([/3]) = 5X](A.)' = l. 



(32) 



fc=i 



Since J([/3]) is continuous and the set given by (13211 is compact, function J([/3]) attains minimum. 



(1) fl(l)^ 



denoted by AL , at some point [Z?*-^^] = (Z?} , ...,l3m). If this procedure is carried out for 



1,2,.. ., we obtain a sequence of numbers A^ , Aj , Because A^:!|_;^ < Xhi and J(]j) is bounded 

from below, we can find the limit 

hm AW=A(i). 



Step 2. Let 



y 



(1) 



{x) = —= 2, Pk sin(A:a::) 



k=l 



denote the linear combination (|30p achieving the minimum AL • We shall prove that sequence 
(j/m )mGN contains a uniformly convergent subsequence. From now on. for simplicity, we will write 
ym instead oiym . Recall that 



p{x){^Dl^y„:f + q{x)yl 



dx 



is convergent, so it must be bounded, i.e., there exists constant A/ > such that 



p{x){:Dl^y„;)+q{x)yl 



Therefore, for all ttt, e N it holds the following: 



dx < M, m e N. 



p{x){'-Dl^y„,)^ dx<M + 



q{x)yl^dx 



< M + max 

a:e[0,7r] 



q{x) 



Waix) 



Wa{x)y^ dx = AI + max 

2:e[0,7r] 



q{x) 



Wa{x) 



Moreover, since p{x) > we have 



Ml. 



mm pix) i'D^^y^) dx < p{x) {^ D^ ^y,„) dx < M^, 
2:e[0,7r] Jo ' Jo 



and hence 



r {^'D^^+y^f dx< 
Jo 



Ml 



min p(x) 



=:M2 



(33) 



Using ([5]), (|33p . condition i/m(0) = and Schwartz inequality, we obtain the following: 

2 



\y.n{xf = |/o" + ^i??.+2/,„(x)|' ^ 



r(«) 



{x^ir-^^Dl + yra{i)dt 



< 



1 



r(a) Jo r(a) 2a - 1 



.2a-l 



so that (ym)meN is uniformly bounded. Now, using Schwartz inequality, equation (|33|1 and the 
fact that the following inequality: 

Vti > t2 > 0, {h - ta)^ < i? - ^2 

holds, we have for any < xi < a:;2 < tt that 

\ymix2) - ym{xi)\ = |/o"_^'=i:»^_^2/™(x2) - II^^+''D^^^y^{xi)\ 
1 



Via) 



(X2 - t)"-^^i?o"+y„,(t)dt - / (xi - t)"-^^i?;^,+y,„(i)dt 



1 



r(«) 



< 



Jo 

(X2 - tf^-'^dt] ' f /■'' [(^^oVy™W) 

'^o,+2/™(t)) 



r(a) 



dt 



+ / {ixi^tr~'-ix2-t)"~Ydt 



dt 



< 



/M2 
Via) 



(^2 _ t)2(«-l)df + / ((xi - t)2("-l) - (X2 - tf^'^-'A dt 



/M2 



r(a)V2a-l 



(X2 - Xi)"-^ + [(X2 - Xl)2"-1 - X^"'! + 



X?"'l+x2"-ll^ 



< 



2 T~ -^1 



T{a)^/2^^~^ 



{X2 - X-i)"- 



Therefore, by Ascoli's theorem, there exists a uniformly convergent subsequence (j/m„)neN of 
sequence (?/m)mGN. It means that we can find y'-^' G C[0,tt] such that 

y^^'> = lim yjn„- 

A'iep 3. Observe that by the Lagrange multiplier rule at [/3] = [/J*-^-'] we have 

d 







df3, 



J([/3]) - AW/([/3]) |[^]=[^a,], j = l,...,m. 



Multiplying each of equations by an arbitrary constant C^ and summing from 1 to m we obtain 



d 



i=i 



d/3, 



= E^'7)^ J([/3])-AW/([/3]) 



l[/3] = [/3< 



(34) 



Introducing 



/im(3:;) = , 7 C-^ sm{jx) 



j=i 



10 



we can rewrite p4p in the form 



= 



pixYD^^yrn^D^^hrn + [^(a;) - \^^Wa{x)]ymK 



dx. 



(35) 



Using the differentiation properties and formula '''D^^y,n = DQ^y,n = DlQ_^_°'y,n we write psp as 



= 



JQ 



+p{x)I^^°'y„,''D^^K„\lZo + / [9(2;) - A(;i)u;a(x)]y„/i„d2: := /„ 



(36) 



By Lemma m (with w ~ l/y/Wa) and Lemma [5] (Appendix) , for function h fulfiUing assumptions 
of Lemma m we shall obtain in the limit (at least for the convergent subsequence (ym„ )neN) the 
relation 

dx+ 



-p'{x)ll-^y^^) '■'D^^h-p{x)I^-"y^^^DWS^h 
+ I [q{x) ~ X'-^^Wc{x)]y^^^hdx -.^ I. 



= 



Let us check the convergence of integrals ([55)1 explicitly 

11,.. - I\ 



(37) 



(38) 



< / I - p'ix)I^+''y,n'D^+h„, +p'(x)/i;"y(i) -D^+h\dx+ 
\p{x)I^-''y„,D-D^+h,n - p{x)I^+"y'^^^ D-D^^h\dx+ 

+ |p(x)/i;"y„,^Z?^+;i„U=o -p(a;)/o;"2/^'^ ^i??+/i|.=o| + 
+ |p(x)/i;"y„,^i5^+/i„,U=^ -p(x)/oi;"2/(i) ^i??+/j|.=.| + 

+ / Mx) - X^^^w^{x)]y^h^ - [q{x) - \'^^^w^{x)]y'^^^ h\dx. 
Jo 

For the first integral we get 

\-p\x)I^^''yrn'-'D^^h„,+p'{x)I^-"y^''> 'D^^hldx < 

< Ib'll • [W'Dl+hW ■ ||/o\r(y™-y''^)IU^ +^'^3i^i-oV^II'^o+(/i™-/i)||L2 

where constant A/3 — sup ||ym|| and || • || denotes the suprcmuni norm in the C[0, tt] space. Now, 

jriGN 

we estimate the second integral 

\p{x)I^+''yr^D-D^_^h,n~p{x)I^+''y^^^ D-D^^h\dx < 



< 

For the next two terms we have 



s we have 



(39) 
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resulting from the convergence of sequence \\ym — y\\ — > 0. For sequence hm ~ gm/y/Wa, we infer 
from Lemma [5] that 

Um \\h'-h'\\^0. 



m ^oo 



Hence, also 



hm \\-D^ + %-,,-h)\\= Um \\ll-+"{K,-h')\\=Q 



and at points x = 0, tt wc obtain 

The above pointwise convergences (|5U|) and (|1(I)) imply that 

Finally, for the last term in estimation ([55)1 we get 

Mx) - \[l^w^{x)]y^h,^ - [q{x) - A(i)u;„(a;)]y(i) h\dx < 



(40) 



lim 

m ^oo 



lim 

m >QO 



TT 



■K 



< I \qix){y^h^ - y(i) h)\dx + / \wa{x){\i\'>y,nh^ - A^^^y^^' h)\dx < 
h Jo 



— "" ' Ikl 



M.-\\h^-h\ 



■ WVm-y 



(1)1 



+n ■ Wwo^W ■ [a (m3 • \\hm - h\\ + \\h\\- ||y,„ - y(i)||) + llyW/ill • |A« - A^l 
where constants M3 = sup \\ym.\\ and A = sup |Am |- We conclude that 

mGN meN 

0= lim /,,„ = / 



ni > 00 



and ([37]) is fulfilled for function y^^^ being the limit of subsequence {yni„) of the sequence {ym)„i£t> 
Step 4- Let us denote in relation ([37|) : 



7i(x):=(Q(.T)-A(i)z«„)yW(x), 

73(x):=-p(a;)/oV"2/(^Hx). 
We observe that jj e C[0, tt], j = 1, 2, 3 and Dq^^j^ e -^^[0, tt] because 

n^^^-vo — n^^" Tt) • t^^'^iMA — t" n (n- 7^i^"i/(i) 
-^^0+ 73 — -^^0+ i^p -'0+ y J — -'0+-^ (^p -'0+ y 

— T" iri ■ 7^1""7;(1' + n • '^ D" 7/(1' 

Both parts of the above function belong to the L^[0, tt] space. 

Assuming that function h in ([37]) is an arbitrary function fulfilling assumptions of Lemma [3] and 
applying Lemma [3] part (1), we conclude that 73 = —p ■ /q^"i/(i' G Ci[0,7r] . From this fact it 
follows that p ■ Dlg^^y^^' <E C[0,7r] and integral ([37|) can be rewritten as 







pix -D^+y^'^'D^+h + (<z(.t) - A(i)«;„(x))y(i)/i' 



dx. 
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Now we apply Lemma |3] part (2) defining 

71 (a;) := 71 (a;) = {q{x) - X^'^w^)y^^\x), 

This time 712 G C'iOjTr] and from Lemma[3]part (2) it follows that 

[-DZMxYD^^+q{x)\ y(i)(x) = A(i)«;„(a;)2/(i)(x). 
By construction this solution fulfills the Diriehlet boundary conditions 

2/(i)(0)=y(i'W=0 
and is nontrivial because 

/(2/^'))= r wcix) (y^'H^))' dx = 1. 



'0 
In addition, we also have for the solution 

Let us observe that from the Diriehlet boundary conditions it follows that y'^^^ also solves the 
FSLP digi-dlTl) in [0,7r]. 

Step 5. Now, let us restore the superscript on y™ and show that [ym ] itself converges to 

V / meN 

y(i). First, let us point out that for given A the solution of 

[^D'^_pixrDS+ + qix)] y{x) = \w^{x)y{x), (41) 

subject to the boundary conditions 

y(0) = y(7r) = (42) 

and the normalization condition 

/ Waix)y'^ dx^l (43) 

Jo 

is unique except for a sign. Next, let us assume that y'^^ solves Sturm-Liouville equation (j41|l 
and that corresponding eigenvalue is A = A'-'^'. In addition, suppose that y^^' is non trivial i.e., we 
can find xq € [0,7r] such that y^^-'(xo) 7^ and choose the sign so that y^^\xo) > 0. Similarly, for 
all m e N, let yin solve (|4T|) with corresponding eigenvalue A = Am and let us choose the signs 
so that yin (xq) > 0. Now, suppose that [ym ] does not converge to y'-'^'. It means that we 

can find another subsequence of ( y™ ) such that it converges to another solution y'^^ of (14ip 

V / mSN 

with A = A^^). We know that for A = A'^^ solution of jH]) subject to (jl2| and ([43]) must be unique 
except for a sign, hence 



yW = _y(l) 



and we must have y'^-'(xo) < 0. However, it is impossible because for all mGN value of yin in xq 
is greater or equal zero. It means that we have contradiction and hence, choosing each yfn with 
adequate sign, we obtain y,„ -^ V^^^ ■ 

Step 6. In order to find eigcnfunction y^^' and the corresponding eigenvalue A*^^' , we again minimize 
functional ([28|) subject to (|29|1 and ([27|l . but now with an extra orthogonality condition 



Wa{x)y{x)y'-^^ (x) dx = 0. (44) 





If we approximate solution by 



ym{x) = —== y^ /3k sm{kx), ym(0) = ym(7r) = 0, 

V " k=l 
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then wc again receive quadratic form (|3ip . However in this case admissible sohitions are points 
satisfying ([5^ together with 



?E/5^4-^=0, (45) 



2 

fc=i 

i.e., they lay in (771— l)-dimensional sphere. As before, we find that function J([/3]) has a minimum 
\m and there exists A^^^ such that 

A(2) = hm A(„2\ 

rn— ^00 

because J{y) is bounded from below. Moreover, it is clear that the following relation: 

A(i) < A(2) (46) 

holds. Now, let us denote by 



£^)(^) = -LX]/3[.'^sin(M, 

the linear combination achieving the minimum Am , where /3'^' ~ {P\ ,■■ ■ ,13m) is the point sat- 
isfying (j32p and (|45p . By the same argument as before, we can prove that the sequence (ym )meN 
converges uniformly to a limit function y"% which satisfies the Strum-Liouvillc equation (|26p with 
A'-^\ the boundary conditions (P7|) . normalization condition (pS)) and the orthogonality condition 
(|44)) . Therefore, solution y'-^^ of the FSLP corresponding to the eigenvalue A'^' exists. Further- 
more, because orthogonal functions cannot be linearly dependent, and since only one eigenfunction 
corresponds to each eigenvalue (except for a constant factor), we have the strict inequality 

Ad) < A(2) 

instead of (|46p . Finally, if we repeat the above procedure, with similar modifications, we can 
obtain eigenvalues A''^', A'"'-', . . . and corresponding eigenfunctions y^^\y^'^\ .... D 

2.2 The First Eigenvalue 

In this section we prove two theorems showing that the first eigenvalue of problem ([26l) - ((27)) is a 
minimum value of certain functionals. As in the proof of Theorem [2] in the sequel, for simplicity, 
we assume that a = and & = tt in the problem (I26 p - ((27|) . 

Theorem 3. Let yd) 5g ^/jg eigenfunction, normalized to satisfy the isoperimetric constraint 

I{y)= I w^{x)y''dx^l, (47) 



associated to the first eigenvalue Ad) of problem (|26p - (|27p and assume that function D"_{p{xYDQ^y) 
is continuous. Then, yd) {g a minimizer of the following variational functional: 

J{y) = f [p{x){'^D^+yf + q{x)y^] dx, (48) 

Jo 

in the class C[0, tt] with '^D'^j^y G C[0,7r] subject to the boundary conditions 

y(0) = y(7r) = (49) 

and an isoperimetric constraint (J47]). Moreover, 

J(yW)-AW. 
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Proof. Suppose that y £ C[0,tt] is a minimizcr of J and '^D^^y e C[0,7r]. Then, by Theorem [I] 
there is number A such that y satisfies equation 

[DZMxfD^+ + qix)] V{x) = \w^{x)y{x), (50) 

and conditions ([47|) . ([49|) . Since D'^_(p(xY W^j^y) and '^D'^_(p{xYDQ^y) are continuous, it foUows 
that p{xY DQ^y{x) \x=tc — 0. Therefore equation ([50|) is equivalent to 

['D'^^xYD^^ + g(a;)] y{x) = \w^{x)y{x). (51) 

Let us multiply (|50p by y and integrate it on the interval [0,7r], then 

{y ■ D'^_{p{xYD^+y) + q{x)y^) dx ^ X f w^ix)y^ dx. 



Applying the integration by the parts formula for fractional derivatives (cf. (|12p ) and having in 
mind that conditions ([i^]) . (|T7)) and p{xYD^+y{x) \x=b = hold, one has 



{-Dl+yfp{x)+q{x)y') dx = X. 



Hence 

Jiy) = X. 

Any solution to problem ([Tf|) - (P^ which satisfies equation ([FT]) must be nontrivial since (|47p 
holds, so A must be an eigenvalue. Moreover, according to Theorem [2] there is the least element in 
the spectrum being eigenvalue A^^^ and the corresponding eigenfunction y^^' normalized to meet 
the isoperimetric condition. Therefore J{y^^'>) = A^^^. D 

Definition 4. We will call functional R defined by 

Jiy) 



Riy) 



livY 



where J{y) is given by (|48p and I[y) by (j47p . the Rayleigh quotient for the fractional Sturm- 
Liouville problem (j26p - (|27p . 

Theorem 4. Let us assume that function y G C[0,7r] with '^D'^j^y G C[0,7r], satisfying boundary 
conditions y{0) ~ J/(7r) ~ and being nontrivial, is a minimizer of Rayleigh quotient R for the 
Sturm-Liouville problem (|26p - (|27p . Moreover assume that function D'^_{j){xYD^j^y) is continu- 
ous. Then, value of R in y is equal to the first eigenvalue X^^' i.e., R{y) = X^^' . 

Proof. Suppose that function y S C[0,7r] with '^D^j^y € C[0,7r], satisfying y{0) = y{7T) = and 
nontrivial is a minimizer of Rayleigh quotient R and that value of i? in y is equal to A, i.e.. 

Consider one-parameter family of curves 

y = y + hij, \h\ < e, 
where r] e C^[0,7r] is such that 77(0) = i]{tt) = 0, tj Y a-nd define the following functions 

0: [-£,e] — > R 



h I — > I{y + hr]) ~ / Wa{x){y + hr]) dx, 
Jo 



ip: [-e,e] 



h ^ J{y + hi^) = / [pix)YD^4y + hTj)Y + q{x){y + hr,f] dx 
Jo 
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and 



C: he,; 



ny+h'i) 



Since C is of class C^ on [— e,e] and 



h ^ R{y + hr,)-^§m 



m<C{h), \h\<e, 



wc deduce that 



Moreover, notice that 



C'(0) = ^Riy + hr,) 



h=0 



<:'<'" -M (*'""- ii*'"" 



and that 



V''(0) = ^J{y + M 



= 2 r [pix) ■ ^D^^y ■ ^D^^rj + q{x)y7j] dx, 

=0 JO 



/''(o)- ^Hy + hv) 



= 2 / [wa{x)yT]] dx. 
h=o Jo 



Therefore 



C'(0) - ^^Riy + hrj) 



h=0 



i{y) 



[p{x) ■ -D^+y ■ -D^^7^ + q{x)yf^] dx ^ ^ J K(a;)OT] dx]^0. 



Having in mind that ^jjH = A and 7/(0) = '/(tt) = 0, using the integration by parts formula ([T! 
we obtain 

/ {[D^_p{xrD^+ + qix)] y{x) - \w^{x)y{x)) r,{x) dx = 0. 

"'0 

Now, applying the fundamental lemma of the calculus of variations we arrive at 

{D^^pixYD'^^ + q[x)\ y{x) = Xw^{x)y{x). (52) 

Under our assumptions p{xY' DQ^y{x) \x=-n = and therefore equation (|5^ is equivalent to 

[^DZM^YD^+ + 9(^)] y{^) = >^wc.{x)y{x). (53) 

Since y 7^ we deduce that number A is an eigenvalue of ([55]) . On the other hand, let A'-™-' be an 
eigenvalue and y^"^^ the corresponding eigenfunction, then 



[^D''^^p{xfD^^+q{x)\ y(")(x) = \'^'^^w^{x)y^"'\x). 
Similarly to the proof of Theorem [31 we can obtain 

J; ({'Dl^y(-^)fp{x)+q{x){y^-^)f) dx 



(54) 



/J'A(™)w«(x)(y('"))2da; 



= A(^ 



for any m G N. That is i?(y'^™^) = '.|^(,„).^ = A'-'"'. Finally, since minimum value of _R at y is equal 
to A, i.e., 



A < -R(2/(")) = A^'") Vm e N, 



we have A = A^^^. 



□ 
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2.3 An Illustrative Example 

Let us consider the following fractional oscillator equation 

[p ■ 'D^_^D^^ - A] yix) = 0, (55) 

where y{a) ~ yih) — and parameter p > 0. One can easily check that problem of finding 
nontrivial solutions to equation ((55|) and corresponding values of parameter A is a particular case 
of problem (|26 p - (|27p with p{x) = p, q{x) = and Wa{x) = 1. The corresponding minimized 
functional is 

J a 

with the isoperimetric condition 

y {x)dx = 1. 



Let us fix the value of parameter p and assume that orders ai,a2 fulfill the condition: ^ < ai < 
Q2 < 1. Then, we obtain for functionals Ja^, Ja^ the following relation 

j„,(y) ^ wvp'D^^^ywi, = wv^ia^'Dywh = \\v^i:r'''iir'Dy\\h 

where we denoted 



T{a2 -ai+1) 
We observe that in the above estimation two cases occur 

Ka2-ai < 1 Jaiiy) < «/a2 (y) 

The relations between functionals for different values of fractional order lead to the set of inequal- 
ities for eigenvalues A^-'-' valid for any j £ N: 

K^2-a,<l A(-'")(ai)<A(^")(a2) 

i^„,_„, > 1 A(-'")(ai) < Kl^^^^ ■ A(^")(a2). 

In particular, when order a2 = 1 we get 

j^Ay) = WVl^'D^iyWl. ^ WV^ilr'DyWh 
< K'!-a, ■ WVpDyWh = Kto^Mv) 

and the following relations dependent on the value of constant A'i_q^ 

Ki^a,<i JaAy)<My) 
i^i_„, >i JaAy)<Ktc.^■Jliy)■ 

Thus comparing the eigenvalues for the fractional and the classical harmonic oscillator equation 
for boundary conditions y{a) ~ y{b) = we conclude that the respective classical eigenvalues are 
higher than the ones resulting from the fractional problem for any j S N, namely 

/^i_„, <1 A«(ai)<A(-'")(l)=p('^') (56) 

X,_>1 AO-)(.0<i^U^AO-)(l)^p( ^^_^^,f;^^,^J ^ (57) 
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3 Appendix 

We shall prove two lemmas, concerning certain convergence properties of fractional and classical 
derivatives, that play an important role in the proof of Theorem[21 Let us begin with the following 
definition of Holder continuous functions. 

Definition 5. Function g is Holder continuous in interval [a, b] with coefficient Q < P < 1 iff 

sup — i i^g — < oo. (58) 

We denote this class of Holder continuous functions as Cfj[a,b]. 

Lemma 4. Let a e (0, 1), functions w,g £ C^[Q, tt] n C\j[—t:, it] he odd functions in [— tt, tt] such 
that w" , g" G L'^[0,it]. If we denote as g^ the m-th sum of the Fourier series of function g, then 
the following convergences are valid in [0, tt] 

lim |ri?o%5™-'^?,+.9llL^-0 (59) 

m —7- oo 

lim \\D'D^+g„, - D'D^^gWL, = (60) 

m - — > oo ' 

lim \\-D^+wg,n-'D^+wg\\L2=0 (61) 

lim \\D'D^+wg,n-D'Dl^wg\\Li=0. (62) 



rn >-oo 



7n > oo 



Proof. We can apply equation (fT3|l and estimate the ||^I?Q_^5^n ~ ^-Dq+^IU^ norm in [0,7r] as 
follows 

\\c r^ct c r^Oi \\ mtI — a// /\\\ ^ t/- \\ > > \\ 

II -Do^+5m - -Do,+5||l2 = lKo,+ (.9m - a)\W < ^1-a • ||g„ - ff ||l2. 

For odd functions from the C^[0,7r] n C]^\—'k,'k] space, g'^ is the yn-th sum of the Fourier series 
of the derivative g' . Hence in interval [0,7r] 



lim ||5:„-5'IIl^-0 

m >-oo 

and from the above inequalities it follows that ([59| is valid in [0, tt] 

lim |ri?o%5,„-^i5o%.g|U2=0. 

Tfl fOO 

Let us observe that for x > 



D^'BZ,^g(x) = DZ^^g\x) = 'Dl^g'ix) + 

mate tl 
T3|) for /3 = 1 — a and p = 1 



r(i-a) 

g'(0) • x-" 



r(l-a) ' 
Therefore we can estimate the distance between D'^Dq _^gm and D'^D^^g in interval [0,7r] using 



\\D''Dl^g,,,-D-Dl^g\\L^< 

< ir^o,+(5m - g')\W + ll(.9m(0) - g'm ■ Y^^^\W = 

= II4T(.C - a")\W + 1.9m (0) - .g'(0)| • II ^f^lU^ < 
< A'i_„ -Wg';^- g"\\L. + |5:„(0) - ^'(0)1 • JT^^^ < 
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< Xi_„ -V^-Wg':^ - g"\W- + |5:„(0) - .g'(0)| 



r(2-a)- 
By assumptions we have in [— 7r,7r] (thence also in [0, tt]) 

hm ||5™-9"I|l^=0 hm |g:„(0) - .9'(0)| == 0. 



771 >QO ra >cxD 



Hence we conclude that (|5(I| is valid. 

The convergence given in (|6ip follows from (|59p . namely 



= ll^',r[(^ff™)'-(«^.9)1lU^< 

< ii^',r^(5™ - 9')\w- + ii^\r (^)' (.9>n - .9)iil^ < 

< ifi_„ • |h(5:„ - 5')IIl^ + Ki-o. ■ II («^)' (5777 - .9)||l^ < 

< A'i_a (llwll • ||g™-g'||L2 + 11 (u;)'|| • 11.9777 -.9 1 1 lO , 

where || • || denotes the supremum norm in the C[0, tt] space. From assumptions of our lemma it 
follows that in [0, tt] 



lim ll5m"5'llL2 = lim ||g„ - gjl^a = 0. 



.Jm"5'l. 

777 ^OO 777 ^OO 



Thus convergence (pT|) is valid. 

To prove convergence (|62p we start by observing that for x > 

D^Dl^wgUx) = Dl^ [wgj [x) = ^'Dl^ {wg„,)' {x) + ^""^"^^^^^^'J^'^ 

D^Dl^wg{x) = Dl^ {wg)' (x) = ^i?„% [wg)' (x) + ^""ff] '^"^"^ ■ 
For the i^-distance between D'^Dq ^wgm and D'^Dq _^_wg in interval [0, tt] we have 

\\D'Dl^wg„, - D-D^^wgllLi < (63) 

< |ri?oV [("'.9,77)' - (wg)'] IU7 + II [(«;ff™)' (0) - (wg)' (0)] • jt^^^IIli = 



< ||/o\r [i^9m)" - (wg)"] lUi + |(«;ff777)' (0) - (wg)' (0)| • 



r(i-a) 



lil 



< 



< Ki^a ■ II (w.9777)" - (wg)" IIli + \{wgr,^)' (0) - (wg)' (0)1 • -— < 

T(2~a) 

1 _^ 

< A'i_Q • Vtt- II {wgm)" ~ (wg)" \\l2 + \{wgm)' (0) - (wg)'(0)| 



r(2-a) 

Because 

{wg,n)" - (wg)" = 

= y^ig'Jn - 5") + 2 H' ■ (g'm - ff') + H" ■ idm - g) 

we have 

II {wgm)" - {wg)" Wl-^ < 

^ W II II" "II lO 11/ \'ll II' 'II ill/ \"ll II II 

< Ikll • II5777-5 I|l2 +2- II (if) II ■ ||g™-5 ||l2 + ||(W) ||l2 ■ 113777 -ff||L2- 

From the assumptions of the lemma it follows that for j = 0, 1, 2 

lim Wg'^} -g^HL-^^- 

m J-oo 
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Hence 

lim \\{wgm)" - (wg)" \\l2 =0. 

m >-oo 

In addition 

lim Uwgr,,)' (0) ~ (wg)' (0)1 = 



m ^oo 



= lim |H'(0)(.9™(0)-5(0)) + ^(0)(.g:„(0)-.9'(0))|< 

m >-cx3 

< lim |(u;)'(0)(5™(0)-.g(0))|+ lim |ti;(0)(5;.(0) - .g'(0))| = 0. 

Taking into account estimation (|63p and the above inequalities we conclude that (|62p is valid. D 

Lemma 5. Let a e (^jl) j /3 < a — ^, function w be positive, even function in [— 7r,7r] and 
w' S C^[— 7r,7r]. Function h' is the derivative of h defined by assumptions of Lemma \M and 
formula (|19p . function g is defined as 

g{x) := h{x)w{x). 

If we denote as g^ them-th sum of the Fourier series of function g, then the following convergences 
are valid in interval [0, tt] 



lim ||5:„-.9'||=0 (64) 

lim g'„M^g'{0) (65) 

lim .g™(7r) = g'in). (66) 



m >-oo 



m >OQ 



rn >cjo 



Proof Definition (|19p in interval [0, tt] implies for derivative /i' 

h'{x) = /o\7(a;) + Ax" + Bx^+"', (67) 

where 7 G C[0,7r] and constants A, _B e R are specified by conditions (fT7|) . p^ in the proof of 
Lemma f. Let us observe that x^'^" € C^[0,7r], function x" is Holder continuous in [0,7r] with 
coefficient /3 < a, thus it can be extended to an odd/even, Holder continuous function in interval 
[— tTjTt]. In addition lQ_^'-f{x) is Holder continuous function in [0, tt] with coefficient /3 < a — ^ 
because: 

|:e-2/|'5 - r{a)V2^^T ' ' " r(a)V2^^T 

and can be extended to an odd/even, Holder continuous function in interval [— tt, tt]. Observe that 
for Holder continuous functions in [— tt, tt] we have the absolute convergence of their Fourier series. 
For function g' we obtain in [0, tt] 

g'{x) — h'{x)w{x) + h{x)w' (x). 

Both terms on the right-hand side are by assumption functions from the C|^[0, 7r]-space and can be 
extended to odd/even functions in the C^[— tt, tt] space. Hence their Fourier series are absolutely 
convergent in [— tt, tt]. Concluding, we have for function g' the convergence in interval [— tt, tt] 

lim ||5:„-5'||==0, 



m ^00 



where || • || denotes the supremum norm in interval [— tt, tt]. Thus the sequence g'^-^ of partial sums 
is also absolutely convergent in interval [0, tt]. Formulas (|65|66p are a straightforward consequence 
of this fact. D 
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